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Layered transition metal trichalcogenides with the chemical formula ABX 3 have attracted recent 
interest as potential candidates for two-dimensional magnets. Using first-principles calculations 
within density functional theory, we investigate the magnetic ground states of monolayers of Mn- 
and Cr-based semiconducting trichalcogenides. We show that the second and third nearest-neighbor 
exchange interactions (J 2 and J 3 ) between magnetic ions, which have been largely overlooked in 
previous theoretical studies, are crucial in determining the magnetic ground state. Specifically, we 
find that monolayer CrSiTe 3 is an antiferromagnet with a zigzag spin texture due to significant con¬ 
tribution from J 3 , whereas CrGeTe 3 is a ferromagnet with a Curie temperature of 106 K. Monolayers 
of Mn-compounds (MnPSs and MnPSes) always show antiferromagnetic Neel order. We identify 
the physical origin of various exchange interactions, and demonstrate that strain can be an effective 
knob for tuning the magnetic properties. Possible magnetic ordering in the bulk is also discussed. 

Our study suggests that ABX 3 can be a promising platform to explore 2D magnetic phenomena. 

PACS numbers: 75.70.Ak,75.50.Pp,75.50.Ee,75.50.Dd 


I. INTRODUCTION 

Recent years have seen a surge of interest in two- 
dimensional (2D) atomic crystals due to their highly tun¬ 
able physical properties and immense potential in scal¬ 
able device applications . 1 5 Since the initial discovery of 
graphene, the family of 2D crystals has grown consider¬ 
ably with new additions such as boron nitride and tran¬ 
sition metal dichalcogenides. The emergence of transi¬ 
tion metal compounds in the landscape of 2D crystals is 
particularly advantageous as it opens the door to many 
physical properties not available in graphene . 6-11 For ex¬ 
ample, in monolayer M 0 S 2 the large spin-orbit interac¬ 
tion leads to a unique spin-valley coupling which might 
be useful for spintronic applications. Another inter¬ 
esting possibility brought by transition metal elements 
is magnetism , 17 a property still missing in the current 
line-up of 2D crystals. In this regard, transition metal 
trichalcogenides (TMTC) such as MnPSs represents a 
rather attractive material family. Similar to dichalco¬ 
genides, these are layered compounds with weak inter¬ 
layer Van der Waals interactions. Furthermore, these 
materials are known to exhibit a large variety of magnetic 
phases , 18-23 making them ideal candidates for exfoliated 
2D magnets. The successful fabrication of a truly 2D 
magnet would also significantly advance our understand¬ 
ing of low-dimensional magnetism. 

Despite the obvious interest and more than three 
decades of experimental studies of bulk TMTC, mag¬ 
netism in these materials remains to be fully understood. 
In particular, even though the spin wave measurement by 
inelastic neutron scattering has pointed out the impor¬ 
tance of exchange interactions beyond nearest-neighbor 
(NN) spins. • relatively little is known about the na¬ 
ture of these interactions and their effect on the magnetic 
ground state. Additionally, the 2 D confinement of elec¬ 
trons upon exfoliation often leads to properties quite dif¬ 


ferent from the bulk crystals. It is thus interesting to ask 
whether there is any change of the magnetic ground state 
when these materials are thinned down to monolayers. 

With these questions in mind, we investigate in this 
paper the magnetic ground states of monolayers of 
Mn- and Cr-based semiconducting TMTC, using first- 
principles calculations within the framework of density 
functional theory (DFT). The Mn-compounds (MnPSs 
and MnPSes) are known to exhibit antiferromagnetic 
(AF) Neel order in their bulk form , 18,25 and are chosen 
here as benchmark for our calculations due to the exten¬ 
sively available experimental data. Interesting proper¬ 
ties such as coupled spin and valley degrees of freedom 
have also been predicted for monolayers of these mate¬ 
rials. 26 The Cr-compounds (CrSiTes and CrGeTes), on 
the other hand, are reported to be ferromagnetic (FM) 
in bulk, thus present a highly interesting system 
to realize 2 D ferromagnets. One of the motivations of 
the present paper is to provide some quantitative under¬ 
standing of magnetism in these compounds. 

Our main findings are summarized below. The major¬ 
ity of the paper is focused on monolayers. We show that 
the second and third NN exchange interactions (J 2 and 
J 3 ) mediated through the p states of chalcogen anions 
are crucial in determining the magnetic ground state. 
Specifically, we find that monolayer CrSiTe 3 is an an¬ 
tiferromagnet with a zigzag spin texture due to signif¬ 
icant contribution from J 3 , whereas CrGeTe 3 a ferro¬ 
magnet with a Curie temperature of 106 K. This result 
is in sharp contrast with previous theoretical studies in 
which only the NN exchange interaction (Ji) was con¬ 
sidered . 30,31 We discuss the physical origin of various ex¬ 
change interactions, and demonstrate that strain can be 
an effective knob for tuning the magnetic properties. A 
uniform in-plane tensile strain of ^ 3% can tune the mag¬ 
netic ground state of CrSiTe 3 from zigzag to ferromagnet, 
with a critical temperature of 111 K. We also find that 
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FIG. 1. (Color online) Crystal and magnetic structure of 
transition metal trichalcogenides ABX 3. The crystal struc¬ 
ture (a) and the top view (b) of monolayers of ABX 3 . The 
transition metal A atoms form a honeycomb structure with 
B 2 Xq ligand occupying the interior of the honeycomb. Top 
view of the different spin configurations: the FM ordered (c), 
AF-Neel ordered (d), AF-zigzag ordered (e), AF-stripy or¬ 
dered (f), with only the transition metal ions shown. Up 
(down) spins are represented by black filled-in (open) circles. 
The crystal structure is drawn using VESTA. 


in bulk CrSiTe3, the intralayer magnetic ordering is very 
sensitive to the out-of-plane lattice constant. For the 
experimental out-of-plane lattice constant (21.0 A), the 
intralayer magnetic ordering is FM in nature. However, 
the interlayer coupling favors AF over FM coupling. This 
is in contradiction with experimental results. 5 One 
possible reason for this discrepancy is discussed, but the 
actual mechanism for ferromagnetism in bulk CrSiTe3 
remains an open question. 


II. CRYSTAL AND MAGNETIC STRUCTURE 


Transition metal trichalcogenides with the chemical 
formula ABX 3 are layered compounds with the struc¬ 
tural space group of R 3 , except MnPS3, which forms 
monoclinic crystals with the C2/m space group. In all 
compounds, the different layers are held together by weak 
van der Waals force. It has been predicted that the 



FIG. 2 . (Color online) The ground state magnetic phase 
diagram for our spin model as a function of J1/J3 and J2/J3. 
Since our calculation finds J3 to be always AF, only J3 > 0 is 
considered. Spins are treated as classical degrees of freedom. 
All compounds studied are located at corresponding param¬ 
eter values. Open symbols are positions under tensile strains 
with arrows indicating the change from the unstrained cases. 


monolayer form of these materials are indeed stable, 30,33 
making them attractive candidates for 2 D magnets. Fig¬ 
ure 1(a) and (b) show the crystal structure of TMTC 
monolayers. The magnetic ions ( A ) form a honeycomb 
lattice within each layer, and each of them is octahedrally 
coordinated by six X atoms from its three neighboring 
(B 2 X 6 ) ligands, with the centers of the hexagons occu¬ 
pied by the B 2 groups. 

Similar to the crystal structure, the magnetic structure 
of bulk TMTC also shows 2D characteristics. It can be 
understood as FM or AF coupled 2 D magnetic layers. 
To describe the 2 D magnetic structure, we consider the 
Heisenberg model on a honeycomb lattice, 

H = Y JiSi J 2 Si A' + Y • % • (!) 

(ij) (m (« ij »> 

where Ji,2,3 are the exchange interactions between NN, 
second NN, and third NN spins. Previous studies have 
shown that it is necessary to include both J 2 and J3 to 
fit the spin wave dispersion from inelastic neutron scat¬ 
tering data. 1 ’ In addition, considering only J\ would 
yield either FM or AF-Neel order, while other magnetic 
ground states have been found experimentally. To com¬ 
pute the exchange interactions, we consider the following 
four possible magnetic ground state: FM, AF-Neel, AF- 
zigzag, and AF-stripy, as shown in Fig. 1 (c)-(f). The 
ground-state phase diagram for our model in Eq. ( 1 ) is 
shown in Fig. 2 as a function of J1/J3 and J2/J3. Here 
J3 is assumed to be positive, as it turns out to be the 
case for all the compounds we studied. It is clear that 
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FIG. 3. (Color online) The partial density of states (PDOS) of the ground states of CrSiTe 3 (AF-zigzag), CrGeTe 3 (FM) and 
MnPS 3 (AF-Neel) are shown in (a), (b) and (c) respectively. These are the three unique ground states exhibited by the ABX 3 
compounds. The PDOS of A, B and X are shown using red lines, broken green lines and dotted blue lines respectively. We 
observe considerable hybridization between the transition metal ( A ) atoms and the chalcogen (A) atoms. 


not only Ji, but also J2 and J3 play an important role 
in deciding the magnetic ground state. 


III. EXCHANGE INTERACTIONS 
A. Computation details 

With the above observation, the magnetic ground 
states of ABX 3 compounds are examined using DFT 
employing the projector augmented wave 34 method 
encoded in Vienna ab initio simulation package 36 with 
the generalized gradient approximation in the parame¬ 
terization of Perdew, Burke and Enzerhof. v We use 
Hubbard U terms (4 eV for Cr and 5 eV for Mn) 39,40 
to account for strong electronic correlations as suggested 
by Dudarev et al Our results were qualitatively insen¬ 
sitive to the different C/’s chosen (2 eV, 4 eV) for the 
Cr-compounds. For each slab a vacuum region more 
than 15 A was used. A cutoff energy of 400 eV and 
a Monkhorst-Pack special k -point mesh of 24 x 14 x 1 for 
the Brillouin zone integration was found to be sufficient 
to obtain the convergence. Structural optimizations were 
performed by fixing the in-plane lattice constants to that 
of the theoretical bulk lattice constants (see Table I). All 
ions were then relaxed with the relaxation of the elec¬ 
tronic degrees of freedom accurate to up to 10 -6 eV. 


B. Exchange interactions 

For each compound, we optimize the crystal structure 
for all four spin configurations [see Fig. 1 (c)-(f)] to 
find the lowest-energy state. Figure 3 shows the partial 
density of states (DOS) of three representative ground 
states: AF-zigzag (CrSiTes), FM (CrGeTes) and AF- 
Neel (MnPSs). It is evident that there is considerable 


hybridization between the chalcogen p states and the 
transition metal d states, further confirming the neces¬ 
sity to include the second and the third NN interaction 
into consideration. By integrating the partial DOS in the 
transition metal atoms for the lowest energy spin config¬ 
uration we obtain Si = 2.45 for Mn-compounds and 2.10 
for Cr-compounds, with the spins having a variation less 
than 0.01 between the different spin configurations. 

To further extract the J’s, we chose to fix the lattice to 
that of the most energetically favorable spin configura¬ 
tion and computed the energies for different spin config¬ 
urations. The exchange coupling constants were derived 
by mapping the DFT energies to the Heisenberg spin 
Hamiltonian ( 1 ), 

^FM/Neel = ^0 + (±3Ji + 6J 2 ± 3J 3 ) |A| 2 , 

-^AF—zigzag/stripy -^0 V (i Jl 2 J 2 ~F ^ J3) \$\ 5 

where E$ is the ground state energy independent of the 
spin configuration. Using these J’s, we also calculated 
the critical temperature by performing a Monte Carlo 
simulation of an Ising model on the 2 D honeycomb lat¬ 
tice. 42 

The magnetic ground state, the computed J’s, along 
with the critical temperature for each compound are 
listed in Table I. The locations of the ground state of 
all the compounds studied are labeled in the phase dia¬ 
gram in Fig. 2. We see that both MnPSs and MnPSe3 
are deep inside the AF-Neel phase. The calculated val¬ 
ues of Ji,2,3 for monolayer MnPSs agree excellently with 
the experimental data for the bulk system, 18 which val¬ 
idates our calculation. We also find that both CrSiTe3 
and CrSiSe3 are in the AF-zigzag phase with the former 
lying close to the boundary of AF-zigzag and FM phase. 
This is different from the FM ground state reported for 
bulk CrSiTes, 27,28 which will be addressed later. Finally, 
of all the compounds studied, CrGeTe3 is the only one 
that has a FM ground state in its unstrained monolayer 
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FIG. 4. (Color online) Top view of monolayer of ABX 3 . Blue, violet and yellow represents A, B and X ions respectively. The 
five possible paths for second NN interaction is show in (a) and the third NN interaction is shown in (b). The two different 
NN hopping paths between two transition metal atoms at different sites (Ai and A 2 ) are shown in (c) with the direct exchange 
as a hopping between the transition metal orbitals and the superexchange interaction characterized by hopping between the 
transition metals through the X atom. 


TABLE I. Lattice constant a, magnetic ground state (GS), 
exchange coupling constants, and magnetic critical tempera¬ 
ture for ABX 3 studied. Critical temperatures are obtained 
from classical Monte Carlo simulations. 



a 

(A) 

GS 

■h 

(meV) ( 

j 2 

meV) ( 

j 3 

meV) 

T c 

(K) 

MnPS 3 (exp) 18a 

5.88 

Neel 

0.77 

0.07 

0.18 

164 

MnPS 3 

5.88 

Neel 

0.79 

0.04 

0.23 

231 

MnPSes 

6.27 

Neel 

0.46 

0.03 

0.19 

147 

MnPSe 3 ( 2 % strain) 

6.40 

Neel 

0.33 

0.03 

0.16 

115 

CrSiSe 3 

6.29 

Zigzag 

-0.74 

0.0 

0.43 

92 

CrSiTe 3 

6.84 

Zigzag 

-1.63 

0.08 

0.71 

160 

CrSiTe 3 (1% strain) 

6.91 

Zigzag 

-1.82 

0.07 

0.66 

130 

CrSiTe 3 (2% strain) 

6.98 

Zigzag 

-1.99 

0.07 

0.60 

72 

CrSiTe 3 (3% strain) 

7.04 

FM 

-2.16 

0.05 

0.54 

111 

CrSiTe 3 (4% strain) 

7.11 

FM 

-2.29 

0.05 

0.50 

158 

CrGeTe 3 

6.91 

FM 

- 1.88 

0.20 

0.22 

106 


a Inelastic neutron scattering was used experimentally to extract 
the energies. 


form. 

We note that J3 is significantly large. This corrob¬ 
orates the decision to include more than just the NN 
interaction. Ignoring it (and J2) had previously yielded 
a different ground state (FM) for CrSiTe3 monolayers 
in previous studies. 30,31 In our calculation for monolayer 
CrSiTe3, FM is indeed lower in energy than AF-Neel. 
The energy difference between the two magnetic states 
was found to be similar to what was reported in Li et 
a/. 30 , when the same U was chosen. But crucially, we 
find that AF-zigzag is even lower in energy than FM, 
and hence is the magnetic ground state. 

Now that we have shown that the interactions have to 
be included up to the third NN not only to interpret the 
neutron diffraction data, but also to get the correct mag¬ 
netic ground state, the imminent task is to understand 
the microscopic origin of the different J’s. We first note 
that J2 and J3 are always AF. Furthermore, the value of 


(a) J 2 : AF 
Top layer 
Middle 
Bottom 



(b) J 3 : AF 
Top 
Middle 



(c) J^AF 



FIG. 5. (Color online) Main contributions of virtual electron 
excitations to magnetic interactions. Virtual electron excita¬ 
tions from two X anions to a pair of (a) second NN and (b) 
third NN A ions, (c) Direct excitations between neighboring 
TM ions (red dots), resulting in AF Ji. (d) Excitations from 
two orthogonal orbitals on a A anion (black dots) to a neigh¬ 
boring pair of A ions. Because of the Hund coupling acting 
in excited states as indicated by a broken circle, this process 
gives rise to FM J\. The actual sign of exchange interactions 
results from the competition between FM and AF contribu¬ 
tions for Ji. Numbers indicate the typical order of the first 
half perturbation processes for each contribution. 


J2 is found to be smaller than J3. Both of these find¬ 
ings are consistent with previous reports on MnPS3 and 
its Fe derivative, FePS3. 18,24 These observations can be 
understood by analyzing the crystal structure. Figure 4 
shows the possible hybridization paths connecting A site 
ions. For the second NN and the third NN A site pairs, 
electrons hop through two X anions [Fig. 4 (a) and (b)]. 
For this reason, J2 and J3 might be regarded as super¬ 
superexchange interactions. Based on the geometry and 
the X anion p states, we expect J2 to be weakly AF be¬ 
cause it involves small X-X hybridizations [Fig. 5 (a)]. 
On the other hand, J3 involves two X anions on the same 
plane, either top layer or bottom layer. Hence, there is 
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strong hybridization of the p states [Fig. 5 (b)], resulting 
in a strongly AF J 3 . 

The NN exchange Ji, on the other hand, shows a large 
variation from compound to compound. The variation is 
so large that it even changes the sign going from the Mn 
compounds to the Cr compounds (see Table I). This be¬ 
havior comes from a unique crystal structure, which nat¬ 
urally gives rise to two competing interactions, i.e., the 
direct exchange and superexchange. The direct exchange 
originates from direct electron hopping between the NN 
A sites [see Figs. 4 (c) and 5 (c)]. For the Mn compounds, 
this exchange is robustly AF as the Mn ions are in the 
half-filled high-spin d 5 state. For the Cr compounds, the 
AF direct exchange is weakened by a FM component as 
Cr ions have partially-filled d shell . 43 The superexchange 
interaction is mediated through the X ions [see Fig. 4 

(c) ]. As the A 1 -X-A 2 angle is close to 90°, this interac¬ 
tion is FM . 44,45 It is important to note that for the su¬ 
perexchange interaction two electrons must excite from 
X anion p states to neighboring A d states [see Fig. 5 

(d) ]. Since the electron excitation energy is large for the 
Mn compounds [Fig. 3(c)], reflecting closed d shell on Mn 
ions, the superexchange is expected to play a minor role 
compared with the direct exchange. On the other hand 
for the Cr compounds, the superexchange could play a 
dominant role. 

To confirm the distinct role of the superexchange mech¬ 
anism’s contribution to Ji, we examine the magnetiza¬ 
tion of a chalcogen ion between two ferromagnetically 
coupled transition metal ions as its magnitude is a good 
indication of the strength of the superexchange interac¬ 
tion. Our DFT calculation for the Cr compounds in the 
FM metastable state showed the total magnetization of 
the chalcogen ions is ^ 0.6 ps per unit cell and hence 
significant. On the other hand, for the Mn compounds, 
the magnetization of the chalcogen ions are an order of 
magnitude smaller, which is consistent with our finding 
of an AF Ji, and the presence of a large electron exci¬ 
tation energy from X p to Mn d. The net result is that 
Ji becomes AF for Mn compounds because of the domi¬ 
nance of AF direct exchange over the FM super exchange, 
while the FM superexchange wins over the AF direct ex¬ 
change making J\ strongly FM for Cr compounds (see 
Table I). Hence, depending on the transition metal ions 
involved, a significant competition is expected between 
FM and AF components, which could lead to a plethora 
of magnetic states. 

This competition is further verified by applying a uni¬ 
form tensile strain on MnPSe 3 and CrSiTe 3 . It is im¬ 
portant to note that both direct exchange and superex¬ 
change do not change sign as we strain the system, but 
the former decreases more rapidly than the latter as the 
atomic distances increase by a tensile strain. As a conse¬ 
quence, | Ji| for MnPSe 3 decreases with strain, where as 
| Ji| for CrSiTe 3 increases with strain (see Table I). This 
confirmed the presence of competing exchange interac¬ 
tions. 

This result immediately suggests the possibility of tun¬ 


ing the magnetic ground state using strain. Here, we con¬ 
sider monolayer CrSiTe 3 as our prototype system and use 
strain as a knob to change the different J’s. CrSiTe 3 is 
an ideal candidate for this study as it lies close to the FM 
and AF-zigzag phase boundary. With a tensile strain, J 2 
and J 3 are both expected to decrease in magnitude as 
the atomic distances are increased. While the effect of 
strain on J\ is subtler, it is expected to increase in mag¬ 
nitude for small strains. Not surprisingly, an application 
of ^ 3% strain leads to a magnetic phase transition with 
ferromagnetism becomes the magnetic ground state (see 
Fig. 2 ). Strain also has a direct impact on the critical 
temperatures. Once the FM ground state is realized, 
the critical temperature T c can be further enhanced with 
strain. As shown in Table I, T c goes up to 158 K for ^ 4% 
strain. With this strong dependence of critical tempera¬ 
ture on the applied strain, it might be even possible to 
engineer room temperature ferromagnetic behavior, for 
large values of strain . 31 


IV. BULK MAGNETIC ORDER 

So far we have only considered the magnetic proper¬ 
ties of monolayer TMTCs. One of the important finding 
is that monolayer of CrSiTe 3 has an AF-zigzag ground 
state, whereas in bulk it is reported to be FM from 
neutron scattering experiments . 27,28 To understand this 
change in magnetic structure when we go from monolayer 
to bulk, we calculated the magnetic ground state of bulk 
CrSiTe 3 . We find that for the experimental out-of-plane 
lattice constant (21.0 A), the intralayer magnetic order¬ 
ing is FM in nature and not AF-zigzag. This switching of 
the intralayer magnetic ordering from AF-zigzag to FM 
as we go from monolayer to bulk is very sensitive to the 
out-of-plane lattice constant of the bulk system. If we 
increase the out-of-plane lattice constant to 22.8 A, the 
intralayer coupling prefers AF-zigzag. This is because 
bulk CrSiTe 3 has ABC stacking, thus an in-plane AF- 
zigzag spin configuration costs more energy compared to 
a FM spin configuration when the interlayer exchange 
interaction becomes strong. 

However, we also find that in bulk, AF interlayer 
coupling is preferred over FM interlayer coupling by 
10.6 meV per Cr atom. Analyzing the interlayer in¬ 
teractions between the Cr atoms, it is dominated by 
super-superexchange. Based on the geometry and be¬ 
cause of the presence of large chalcogen atoms, we ex¬ 
pect this interaction to be comparable to J 3 and hence 
significant. As the mechanism for interlayer super¬ 
superexchange coupling is similar to that of the intralayer 
super-superexchange coupling previously discussed for J 2 
and J 3 , it is not surprising that these interactions are AF 
in nature. But this contradicts experimental findings of 
bulk ferromagnetism in CrSiTe 3 . 2/,2S 

One possible source for this discrepancy is the ab¬ 
sence of the dipole-dipole interaction in DFT calcula¬ 
tions, as discussed previously. 47,48 By introducing the 
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spin-orbit coupling, we have confirmed that there is an 
out-of-plane easy axis anisotropy, in accordance with ex¬ 
periments . 2 5 5 With this easy axis, the interlayer FM ar¬ 
rangement can become energetically more favorable than 
the AF arrangement, because of the dipole-dipole inter¬ 
action. 48 We also note that for MnPSs the dipole-dipole 
interaction is negligible due to the AF ordering within 
each plane. This could explain why our DFT results 
show excellent agreement with the bulk experimental re¬ 
sults for the Mn compounds. Nonetheless, the actual 
mechanism for ferromagnetism in bulk CrSiTe 3 remains 
an open question. 

V. SUMMARY 

In this paper, we studied the magnetic properties of 
monolayers of van der Waals transition-metal trichalco- 
genides ABX 3 using density functional theory. In or¬ 
der to understand the rich magnetic behavior observed 
in these systems, we derived local spin models using the 
DFT energy of the magnetic ground state and metastable 
excited states. Because of the extended nature of the 
p state of the chalcogen atoms, second nearest-neighbor 
and third nearest-neighbor interactions are found to play 
significant roles. Specifically, we find that monolayer 
CrSiTe 3 is an antiferromagnet with a zigzag spin tex¬ 
ture due to significant contribution from J 3 , whereas 


CrGeTe 3 is a ferromagnet with a Curie temperature of 
106 K. Detailed analyses on the magnetic interactions led 
us to predict that monolayers CrSiTe 3 can be made ferro¬ 
magnetic with the application of a moderate uniform in¬ 
plane tensile strain of 3%, which is experimentally feasi¬ 
ble. Our studies demonstrate transition-metal trichalco- 
genides ABX 3 are possible candidates for spintronic ap¬ 
plications; especially CrGeTe 3 and strained CrSiTe 3 are 
promising for two-dimensional ferromagnetic semicon¬ 
ductors. The magnetic ordering of bulk CrSiTe 3 , how¬ 
ever, remains an open question. 
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